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NLO QCD corrections to hadronic Higgs production with heavy quarks 
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The production of a Higgs boson in association with a pair of ti or bb quarks plays a very important role at 
both the Tevatron and the Large Hadron Collider. The theoretical prediction of the corresponding cross sections 
has been improved by including the complete next-to-leading order QCD corrections. After a brief description of 
the most relevant technical aspects of the calculation, we review the results obtained for both the Tevatron and 
the Large Hadron Collider. 



1. INTRODUCTION 

The existence of a relatively light Higgs boson 
is both suggested by precision fits of the Standard 
Model (SM) and theoretically required by the 
Minimal Supersymmetric extension of the Stan- 
dard Model (MSSM). Searches at both the Teva- 
tron and the Large Hadron Collider (LHC) will 
play a crucial role in testing this hypothesis and 
in discriminating between different models that 
imply the existence of one or more Higgs bosons. 
In this context, the production of a Higgs boson 
in association with a heavy quark and antiquark 
pair, both ti and bb, plays a very important role. 

The associated production of a Higgs boson 
with a pair of tt quarks has a very distinctive 
signature, and can give the only handle on a di- 
rect measurement of the top quark Yukawa cou- 
pling, perhaps the most crucial coupling in ex- 
ploring the origin of fermion masses. Observing 
pp tth at the Tevatron (^s = 2 TeV) will 
require very high luminosity [1] and will proba- 
bly be beyond the machine capabilities. On the 
other hand, if Mh < 130 GeV, pp — > tih is an 
important discovery channel for a SM-like Higgs 
boson at the LHC (^s = 14 TeV) [2-4]. Given 
the statistics expected at the LHC, pp tth, 
with h bb,T^T~' ,W^W~ will also be in- 
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strumental to the determination of the couplings 
of a discovered Higgs boson [5-9] . Several analy- 
ses show that precisions of the order of 10-15% on 
the measurement of the top quark Yukawa cou- 
pling can be obtained with integrated luminosi- 
ties of 100 fb~^ per detector. Moreover, the com- 
bined measurements of pp — > tih with h bb and 
h — > r+T~ could provide the only model inde- 
pendent determination of the ratio of the bottom 
quark to the r lepton Yukawa couplings [7] . 

The associated production of a Higgs boson 
with a pair of bb quarks has a very small cross sec- 
tion in the SM, and can therefore be used to test 
the hypothesis of enhanced bottom quark Yukawa 
couplings which is common to many extensions 
of the SM, such as the MSSM for large values 
of tan/?. Both the Tevatron and the LHC will 
be able to search for evidence of an enhanced 
bbh production, looking for a final state con- 
taining no bottom quarks (inclusive production), 
one bottom quark (semi-inclusive production) or 
two bottom quarks (exclusive production). An 
anomalously large inclusive Higgs boson produc- 
tion will be a clear signal of new physics. In the 
MSSM with large tan/3 this can be mainly as- 
cribed to an enhanced bottom quark Yukawa cou- 
pling. In general, however, it cannot be uniquely 
interpreted, since other production channels, like 
the leading gg h gluon fusion , can contribute 
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as well. Detecting one or two bottom quarks in 
the final state is enough to remove this ambi- 
guity, which is the source of the interest in the 
semi-inclusive and exclusive production channels, 
in spite of their smaller cross section. The ex- 
clusive measurement corresponds to the smallest 
cross section, but it also has a very reduced back- 
ground and this is the case we will consider in the 
following. The final states can be further catego- 
rized according to the decay of the Higgs boson. 
Existing studies have considered mostly the dom- 
inant Higgs decay channel, h ^ bb [10,2], but also 
li+H- [11,12] and h t+t' [13]. 

In view of their phenomenological relevance, a 
lot of effort has been recently invested in improv- 
ing the stability of the theoretical predictions for 
the hadronic total cross sections for pp, pp — > tih 
and pp,pp bbh. Since the tree level or Lead- 
ing Order (LO) cross section is affected by a very 
large renormalization and factorization scale de- 
pendence, the first order or Next-to-Lcading Or- 
der (NLO) QCD corrections have been calcu- 
lated for the inclusive pp,pp tih cross sec- 
tion [14-19] and for the the inclusive and exclu- 
sive pp,pp bbh cross sections [20-24], while 
the inclusive bb ^ h cross section has also re- 
cently been calculated including the Next-to-Next 
to Leading Order (NNLO) QCD corrections [25]. 
In all cases, the NLO cross section has a dras- 
tically reduced renormalization and factorization 
scale dependence, of the order of 15-20% as op- 
posed to the initial 100% uncertainty of the LO 
cross section, and leads to increased confidence 
in predictions based on these results. In this pro- 
ceeding we will present the results of our calcula- 
tion of the NLO cross section for both the inclu- 
sive pp,pp tth [15,16,18,19] and the exclusive 
pp, pp — > bbh cross sections [24] , where h denotes 
the SM Higgs boson and, in the case of bbh, also 
the scalar Higgs bosons of the MSSM. 

The calculation of the NLO corrections to the 
hadronic processes pp, pp tih and pp, pp bbh 
presents challenging technical difficulties, ranging 
from virtual pentagon diagrams with several mas- 
sive internal and external particles to real gluon 
and quark emission in the presence of infrared sin- 
gularities. A general overview of the techniques 
developed and employed in our calculation are 



presented in Section 2, and the corresponding re- 
sults are illustrated in Section 3 for pp, pp tih 
and in Section 4 for pp, pp — > bbh. 

2. THE CALCULATION 

The total cross section for pp QQh (for Q = 
t, b) at 0{al) can be written as: 

(yNLoipP QQh) = ^ ^ I dxxdx2 (1) 

X [^f (xi, (a;2, li)a'Lo{t^) + (1 ^ 2)] , 

where J^f^^^ are the NLO parton distribution 
functions (PDFs) for parton i in a (anti)proton, 
defined at a generic factorization scale /ij = /i, and 
(Tn^o is the C(a^) parton-level total cross section 
for incoming partons i and j, made of the chan- 
nels qq,gg QQh and {q,q)g QQh + {q,q), 
and renormalized at an arbitrary scale /x^ which 
we also take to be Hr = A*- We note that the effect 
of varying the renormalization and factorization 
scales independently has been investigated and 
found to be negligible. The partonic center of 
mass energy squared, s, is given in terms of the 
hadronic center of mass energy squared, Sh, by 

S = XiX^Sh- 

The NLO parton-level total cross section, 
o'nlo^ consists of the 0{al) Born cross section, 
(Tio, and the 0{as) corrections to the Born cross 
section, S&]^,i^o, including the effects of mass fac- 
torization. da^^Lo contains virtual and real cor- 
rections to the parton-level QQh production pro- 
cesses, qq QQh and gg QQh, and the 
tree-level {q, q)g initiated processes, {q, q)g 
QQh{q, q), which are of the same order in Ug. 

The 0{as) virtual and real corrections to qq 
QQh and gg QQh have been discussed in de- 
tail in Refs. [16,19] and we will highlight in the 
following only the most challenging tasks. 

2.1. Virtual corrections 

The calculation of the 0{as) virtual corrections 
to qq, gg QQh (for Q = t,b) proceeds by re- 
ducing each virtual diagram to a linear combi- 
nation of tensor and scalar integrals, which may 
contain both ultraviolet (UV) and infrared (IR) 
divergences. Tensor integrals are further reduced 
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in terms of scalar integrals [26]. The finite scalar 
integrals are evaluated by using the method de- 
scribed in Ref. [27] and cross checked with the FF 
package [28]. The scalar integrals that exhibit UV 
and/or IR divergences are calculated analytically. 
Both the UV and IR divergences are extracted by 
using dimensional regularization in d = 4 — 2e di- 
mensions. The UV divergences are then removed 
by introducing a suitable set of counterterms, as 
described in detail in Refs. [16,19,24]. The re- 
maining IR divergences are cancelled by the anal- 
ogous singularities in the soft and collinear part 
of the real gluon emission cross section. 

The most difficult integrals arise from the IR- 
divergent pentagon diagrams with several mas- 
sive particles. The pentagon scalar and ten- 
sor Feynman integrals originating from these di- 
agrams present either analytical (scalar) or nu- 
merical (tensor) challenges. We have calculated 
the pentagon scalar integrals as linear combina- 
tions of scalar box integrals using the method of 
Ref. [29,30], and cross checked them using the 
techniques of Ref. [27]. Pentagon tensor integrals 
can give rise to numerical instabilities due to tiic 
dependence on inverse powers of the Gram deter- 
minant (GD), GD=det(pi-Pj ) forp, and pj exter- 
nal momenta. wliic;li vanishes at the boundaries 
of phase space when two momenta become de- 
generate. These are spurious divergences, which 
cause serious numerical difficulties. To overcome 
this problem we have calculated and cross checked 
the pentagon tensor integrals in two ways: nu- 
merically, by isolating the numerical instabilities 
and extrapolating from the numerically safe to 
the numerically unsafe region using various tech- 
niques; and analytically, by reducing them to a 
numerically stable form. 

2.2. Real corrections 

In computing the 0{as) real corrections to 
QQ, 99 — > QQh and {q, q)g QQh + (q, q) (for 
Q = t,b) it is crucial to isolate the IR divergent re- 
gions of phase space and extract the correspond- 
ing singularities analytically. We achieve this by 
using the phase space slicing (PSS) method, in 
both the double [31] and single [32-34] cutoff ap- 
proaches. In both approaches the IR region of 
the QQh+g phase space where the emitted gluon 



cannot be resolved is defined as the region where 
the gluon kinematic invariants: 

Sig = ■ Pg = 2EiEg{l " ^ COS 6 ig) (2) 

become small. Here pi is the momentum of an 

external (anti)quark or gluon (with energy Ei), 
Pi = -^/l — mj/Ef, Pg is the momentum of the ra- 
diated final state gluon (quark/antiquark) (with 
energy Eg), and 9ig is the angle between pi and 
Pg. In the IR region the cross section is cal- 
culated analytically and the resulting IR diver- 
gences, both soft and collinear, are cancelled, af- 
ter mass factorization, against the corresponding 
divergences from the 0{as) virtual corrections. 

The single cutoff PSS technique defines the IR 
region as that where 

Sig ^ ^min •> (^) 

for an arbitrarily small cutoff Smin- The two cut- 
off PSS method introduces two arbitrary param- 
eters, 5s and 5c, to separately define the IR soft 
and IR collinear regions according to: 

5s\/s 

Eg < — - — soft region , 

(1 — cos 9ig) < 5c collinear region . (4) 

In both methods, the real contribution to the 
NLO cross section is computed analytically below 
the cutoffs and numerically above the cutoffs, and 
the final result is independent of these arbitrary 
parameters. With this respect, it is crucial to 
study the behavior of a^Lo in a region where the 
cutoff (s) are small enough to justify the analytical 
calculations of the IR divergent contributions to 
the real cross section, but not so small as to cause 
numerical instabilities. 

3. RESULTS FOR tih PRODUCTION 

The impact of NLO QCD corrections on the 
tree level cross section for pp tih (LHC) pro- 
duction in the SM is illustrated in Figs. 1 and 
2. Similar results for the case of pp —t tth (Teva- 
tron) can be found in Ref. [15,16]. Results for cTio 
are obtained using the 1-loop evolution of as{ij) 
and CTEQ5L parton distribution functions [35], 
while results for (Jm^o are obtained using the 2- 
loop evolution of {p) and CTEQ5M parton dis- 
tribution functions, with af^°(Mz) = 0.118. The 
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Figure 1. Dependence of (Jlo,nlo{pP tth) 
on the renormalization/factorization scale ji, at 
/¥;7=14 TeV, for A4 = 120 GeV. 



top quark mass is renormalized in the OS scheme 
and its pole mass is fixed at mt — 17 A GeV. 

Fig. 1 illustrates the renor- 
malization/factorization scale dependence of a^o 
and (Tjvio at the LHC. The NLO cross section 
shows a drastic reduction of the scale dependence 
with respect to the lowest order prediction. Fig. 2 
complements this information by illustrating the 
dependence of the LO and NLO cross sections 
on the Higgs boson mass at the LHC. For scales 
/i > 0.4/io (mo ~ Tit + A//1/2) the NLO correc- 
tions enhance the cross section. We estimate the 
remaining theoretical uncertainty on the NLO re- 
sult to be of the order of 15-20%, due to the left 
over /i-dependence, the error from the PDFs, and 
the the error on the top quark pole mass mj . 

4. RESULTS FOR bbh PRODUCTION 

We evaluate the fully exclusive cross section 
for bbh production by requiring that the trans- 
verse momentum of both final state bottom and 
anti-bottom quarks be larger than 20 GeV {p^ > 
20 GeV), and that their pseudorapidity satisfy the 
condition |ryb| < 2 for the Tevatron and |?7b| < 2.5 
for the LHC. This corresponds to an experiment 
measuring the Higgs decay products along with 
two high pt bottom quark jets. In order to bet- 
ter simulate the detector response, the final state 




100 



120 



140 160 
M, (GeV) 



180 



200 



Figure 2. cr„io(pp 
as functions of Mh-, 
mt -\- Mh/1 and /i = 



tth) and <Jlo{pP tth) 
at y/s^ = 14 TeV, for n = 
2mt + Mh. 



gluon and the bottom/anti-bottom quarks are 
treated as distinct particles only if the separa- 
tion in the aziniuthal angle-pseudorapidity plane 
is AR > 0.4. For smaller values of AR, the four 
momentum vectors of the two particles are com- 
bined into an effective bottom/anti-bottom quark 
momentum four-vector. 

As for tth production, our numerical results for 
the NLO (LO) cross sections are obtained using 
CTEQ5M (CTEQ5L) PDFs [35] and the 2-loop 
(1-loop) evolution of as(/i) . In the bbh case, how- 
ever, we have also investigated the dependence of 
the NLO result on the choice of the renormaliza- 
tion scheme for the bottom quark Yukawa cou- 
pling. The strong scale dependence of the MS 
bottom quark mass {mb{fi)) plays a special role in 
the perturbative evaluation of the bbh production 
cross section since it enters in the overall bottom 
quark Yukawa coupling. The same is not true for 
tth production since the MS top quark mass has 
only a very mild scale dependence. The bottom 
quark pole mass is taken to be mt — 4.6 GeV. 
In the OS scheme the bottom quark Yukawa cou- 
pling is calculated as gf,^;^ = riih/v, while in the 
MS scheme as Shthit'') = ^b{f')/v, where we use 
the 2-loop (1-loop) MS bottom quark mass for 
the NLO (LO) cross section respectively. 

The impact of NLO QCD corrections on the 
tree level cross section for bbh exclusive produc- 
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Figure 3. ctnlo and for pp — > bbh at 

= 2 TeV (top) and for pp — > bbh at ^/s — 
14 TeV (bottom) as a function of the renormal- 
ization/factorization scale ^, for Mh = 120 GeV. 

tion in the SM is summarized in Fig. 3 for both 
the Tevatron and the LHC. In both the OS and 
the M S schemes the stabihty of the cross section 
is greatly improved at NLO, and the correspon- 
dent theoretical uncertainty reduced to 15-20%. 
The MS results seem to have overall a better per- 
turbative behavior, although the variation of the 
NLO cross section about its point of least sensi- 
tivity to the renormalization/factorization scale 
is almost the same when one uses the OS or AIS 
schemes for the bottom Yukawa coupling. This 
indicates that the running of the Yukawa coupling 
is not the only important factor to determine the 
overall perturbative stability of the cross section. 



The difference between the OS and MS results 
at their plateau values should probably be inter- 
preted as an additional theoretical uncertainty. 

Finally, in Fig. 4 we illustrate the dependence 
of the exclusive cross section, at the Tevatron and 
at the LHC, on the Higgs boson mass, both in the 
SM and in some scenarios of the MSSM, corre- 
sponding to tan (3 = 10, 20 and 40. For the Teva- 
tron we consider the case of the light MSSM scalar 
Higgs boson(/i°) while for the LHC we consider 
the case of the heavy MSSM scalar Higgs boson 
{H'^)- We see that the rate for bbh production 
can be significantly enhanced in a supersymmet- 
ric model with large values of tan/3, and makes 
bbh a very important mode for discovery of new 
physics at both the Tevatron and the LHC. 
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